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O Abstract: 

Earlier, two of us and M. Unsal [1] showed that a class of 4d gauge theories, when compactified 
on a small spatial circle of size L and considered at temperatures (3 near the deconfinement 
transition, are dual to 2d "affine" XY-spin models. We exploit this duality to study the 
deconfinement phase transition in SU(3)/1i^ gauge theories with nj > 1 massless adjoint Weyl 
fermions, QCD(adj) on K 2 x §J x S[. The dual "affine" XY-model describes two "spins" - 
Q-i compact scalars taking values in the SU (3) root lattice. The spins couple via nearest-neighbor 

interactions and are subject to an "external field" perturbation preserving the topological Zg 

1—1 and a discrete Zg X subgroup of the anomaly-free chiral symmetry of the 4d gauge theory. 

The equivalent Coulomb gas representation of the theory exhibits electric-magnetic duality, 
which is also a high-/low-temperature duality. A renormalization group analysis suggests — 

^ but is not convincing, due to the onset of strong coupling — that the self-dual point is a fixed 

point, implying a continuous deconfinement transition. Here, we study the nature of the 
transition via Monte Carlo simulations. The Z 3 x Zg X order parameter, its susceptibility, 
the vortex density, the energy per spin, and the specific heat are measured over a range of 

£SJ volumes, temperatures, and "external field" strengths (in the gauge theory, these correspond 

to magnetic bion fugacities). The finite-size scaling of the susceptibility and specific heat we 
find is characteristic of a first-order transition. Furthermore, for sufficiently large but still 
smaller than unity bion fugacity (as can be achieved upon an increase of the §^ size), at 
the critical temperature we find two distinct peaks of the energy probability distribution, 
indicative of a first-order transition, as has been seen in earlier simulations of the full 4d 
QCD(adj) theory. We end with discussions of the global phase diagram in the f3-L plane for 
different numbers of flavors. 
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1. Introduction 

1.1 Motivation: compactification on as a theoretical laboratory 

The observation that the dynamics of certain four-dimensional (4d) gauge theories consider- 
ably simplifies upon compactification on a spatial circle of small size L — i.e., upon considering 
the theories on M 1 ' 2 x instead of M 1 ' 3 — originated in studies of 4d supersymmetric theories, 
compactified on from the late 1990's. Many results, some of which were new, on super- 
symmetric gauge dynamics were obtained by small- L studies, combined with the "power 
of holomorphy" [2]. In 2007, Unsal realized [3] that this simplification transcends super- 
symmetry and also occurs in a large class of nonsupersymmetric 4d gauge theories, notably 
the ones with adjoint massless fermions, QCD(adj). At small L, these theories dynami- 
cally abelianize and many nonperturbative properties, usually not amenable to an analytic 
treatment — confinement, the generation of mass gap, and discrete or abelian chiral symmetry 
breaking — become tractable within a theoretically controlled semiclassical approximation. It 
was found in [3] that in QCD(adj), the nonperturbative mass gap is due to a new kind of 
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topological "molecules" — the "magnetic bions", carrying zero topological but nonzero mag- 
netic charge. The magnetic bions populate the vacuum and cause confinement of electric 
charges by the 3d Debye screening mechanism of Polyakov [4] . It is important to stress that 
while the long-distance confining dynamics is, indeed, three-dimensional, the mere existence 
of the magnetic bions is due to the locally 4d nature of the theory (in the 3d L = limit there 
are no magnetic bions). These small-L studies showed that a piece of conventional wisdom, 
that "confinement is a pure-glue phenomenon," does not hold universally — the nature of the 
topological excitations causing the mass gap crucially depends on the massless fermion con- 
tent. This dependence was later found to be a generic phenomenon, occurring in a variety of 
nonsupersymmetric gauge theories with massless fermions, both chiral and vectorlike [5,6]. 

After the small-L physics of confinement is semiclassically understood, it is natural to 
also study the thermal behavior of the theory. Not surprisingly, one finds that it also has 
interesting and novel features. It is well-known that 4d gauge theories undergo a confinement- 
deconfinement phase transition at some critical temperature. This transition persists upon 
compactification to small-L (the theories are now considered on I 2 x §^ x §[, where (3 is the 
inverse temperature). The small-L finite-T dynamics can also be described semiclassically: 
near the deconfinement transition, the partition function of the gauge theory is that of a 
2d Coulomb gas of electrically and magnetically charged particles, interacting via long-range 
Coulomb (and dual-Coulomb) interactions and by Aharonov-Bohm phase interactions. This 
is similar to the case of the 3d Polyakov model or deformed pure Yang- Mills theory [7], but 
because of the "molecular" nature of the topological excitations causing confinement, the 
dynamics has a different, rather rich, and yet unexplored structure. 

In what follows, we will describe the map between the physics of thermal 4d QCD(adj) 
on and the spin model that we study in this paper (or the equivalent 1 Coulomb gas). We 
will not give a derivation, for which we refer the reader to [1], but provide this discussion in 
order to elucidate the relevant scales and the relation of the quantities of interest in QCD(adj) 
to spin-model observables. We emphasize that the dual spin model is more than a Landau- 
Ginsburg effective theory of an order parameter: owing to the small-L calculability, we have a 
detailed understanding of how the microscopic physics of the gauge theory is reflected in the 
affine spin model. We will also outline the qualitative picture of the deconfinement transition 
as occurring due to a competition between electric and magnetic degrees of freedom. The goal 
of this paper is to study this picture in more quantitative detail, with emphasis on finding 
the order of the deconfinement transition by numerical simulation of the dual spin model. 

1.2 Digression: on the small-L/large-L relation 

Before continuing with the study outlined above, let us comment on the relation between 
small- and large-L physics. The magnetic bion mechanism of confinement in QCD(adj) is 
operative at small but nonzero L <C Aq^ d , where it is under full theoretical control. It 
gives us a remarkable window of calculability permitting the analytical study of difficult 

lr The duality between 2d electric-magnetic Coulomb gases and theories of coupled XY-spins has a long 
history in the condensed matter literature, see [8]. However, the "affine" XY-models dual to QCD(adj), see 
eq. (1.4), especially the ones for N c > 3, are new. 
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nonperturbative phenomena in locally 4d gauge theories — an opportunity, which (we feel) is 
interesting to explore on its own. 

As L increases at fixed A, strong coupling sets in, making it difficult (in theories with- 
out supersymmetry) to continue to the physically interesting case of large L while retaining 
theoretical control. 2 In Seiberg-Witten theory, however, a continuous connection to the 4d 
mechanism of confinement by monopole condensation exhibited by softly-broken N = 2 su- 
persymmetric theories can be made [9]. Thus, an explicit relation exists between the two 
known theoretically controlled mechanisms of confinement in locally 4d nonabelian gauge 
theories (admittedly both mechanisms are essentially abelian). While there should be a con- 
nection to the 4d mechanisms of nonabelian confinement and chiral symmetry breaking (also 
nonabelian) in non-supersymmetric theories as well, we feel that it is not likely to be found 
by purely analytical means (we base this expectation on the complex relation between the 
two mechanisms found in Seiberg-Witten theory and the crucial role played by supersym- 
metry in establishing it). We finally mention two related recent observations. First, the 
use of resurgence and trans-series in quantum field theory [10] may eventually offer a more 
explicit connection between small- and large-L physics, at least in particular theories. Sec- 
ond, it was recently realized that the novel topological "molecules" found at small-L (e.g., 
the already-mentioned "magnetic bions" , the "center-stabilizing bions" , and various kinds of 
monopole-instantons) may be relevant to 4d physics in another guise, namely play a role in 
the microscopic description of the thermal deconfinement transition in nonsupersymmetric 
4d pure Yang-Mills theory [11]. 

1.3 Mapping a 4d gauge theory on to an "afRne" XY model: the melting of 
2d-crystals and SU(3) QCD(adj) 

In this paper, our interest is SU(3)/Zs QCD(adj) with nj massless adjoint Weyl fermions 
on R 1 ' 2 x Ei\. As shown in [1], the theory near the deconfinement transition is described 
by an "affine" XY-spin model. The reason this is possible is that the partition functions of 
both the spin model and thermal QCD(adj) at small-L, for a range of temperatures near the 
deconfinement transition, can be given a 2d Coulomb gas representation. The parameters of 
the two Coulomb gases can then be mapped to each other. In the rest of this Section, we first 
briefly review the continuum dynamics of QCD(adj) and then present the dual spin model 
description. 

1.3.1 Review of the continuum dynamics of QCD(adj) on R 2 x x 

We begin with the zero-temperature case. When the theory 3 is considered at small L, SU (3) 
QCD(adj) dynamically abelianizes. In other words, SU(3>) breaks to U(l) 2 at a scale ~ £ 3> 

2 An exceptional nonsupersymmetric example is provided by theories which at L — > oo run into a Banks- 
Zaks-like weakly coupled fixed point — as the n/ = 5 QCD(adj) is believed to — see [6] and Section 3.3. 

3 In this Section, we will not distinguish between SU(3) and SU(3)/Z3 gauge theory. The distinction will 
become relevant in Section 1.3.2 where we consider the dual spin model, which has the symmetries of the 
SU(3)/2i3 theory. In particular, there is no center symmetry, but a topological Z,\ symmetry (also called "dual 
center symmetry" [12]), see the discussion in Section 1.3.2. 
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h-QCD- Thus, the only perturbative excitations relevant to the dynamics at length scales S> L 
are the two massless photons in the unbroken Cartan subalgebra of SU(3). These can be 
dualized to compact scalars (recall that a photon in 3d is dual to a scalar). Some fermions 
also remain massless (the components of the adjoint fermions along the Cartan generators). 

Nonperturbatively, however, the dynamics is quite complex and various topological ex- 
citations, the magnetic monopole-instantons [4] and twisted [13] monopole-instantons, play 
a central role. The small- L nonperturbative dynamics can be described semiclassically be- 
cause the two U(l) couplings are free and the theory is weakly coupled. The upshot of the 
studies of [3] is that "molecular" (i.e., correlated) instanton events, the magnetic bions, pro- 
liferate in the vacuum and cause the appearance of mass gap and the confinement of electric 
charges (see also [12,14] for recent detailed studies). The nonperturbative mass gap is, up to 
pre-exponential factors: 

^ _ 4tt 2 

m a ~ - e ^> . (1.1) 
Li 

Here, L is the size of §>^ and g\(L) — the 4d gauge coupling, frozen at the scale L. The 
partition function of the zero-temperature theory can be described as a dilute 3d Coulomb 
gas of magnetic bions, interacting via mutual - "magnetostatic" potentials. 

Next, we consider the theory at finite temperature, specifically in the regime: 

e 3 ^m 2 «T<C ^ 4V '- . (1.2) 
L ttL 

The physical meaning of the upper bound on T in (1.2) is that the temperature has to be 
smaller than the inverse size of the bions, so that bions can be treated as pointlike and 
non-dissociated objects, while the lower bound implies that the temperature is larger than 
the inverse of the typical distance between bions, so that the bions are described as an 
effectively 2d gas [1]. It turns out that (1.2) is the interesting range of T for the deconfinement 
transition. At finite temperature, the heavy VF-bosons, of mass ~ ^, resulting from the 
SU(3) — > U(l) 2 breaking, can be excited with probability given by the Boltzmann factor, 

m w O(l) 

~ e t r^j g lt j and their effect on the dynamics cannot be neglected. Since in the 
range (1.2) mw ~ r 3> T, the VF-bosons are essentially static and form a non-relativistic 
2d gas of electrically charged particles. Thus, one arrives [1] at the representation of the 
partition function of the thermal gauge theory as a 2d grand canonical neutral Coulomb gas 
of magnetically and electrically charged particles — the magnetic bions and the VF-bosons. 
The "particles" 4 in each — magnetic or electric — component of the gas have their respective 
mutual (dual) Coulomb interactions, while the "particles" belonging to different components 
have Aharonov-Bohm phase interactions with each other. This two-component Coulomb gas 
of electric and magnetic charges is equivalent to a generalization of the sine-Gordon model 
(see Section 3). 

A qualitative picture of the phase transition is as follows. At low-T, very few VF-bosons 
are present in the magnetic bion plasma, so the confinement of electric charges in the magnetic- 
bion plasma persists, albeit with some amount of screening of the string tension. In this, 

4 To explain the use quotation marks: the electric objects, the W-bosons, are genuine particles, while the 
magnetic bions are instantons. Both appear as particles in the finite-T dimensionally reduced description. 
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T < T c , phase the few VF-bosons present are confined in neutral dipole pairs. As T increases, 
more dipole pairs appear, causing further screening of the confining string tension and thus an 
increase of the size of the dipoles, until, at T c , a transition to a deconfined phase occurs — one 
where W^-bosons form an electric plasma and the magnetic bions are confined in magnetically 
neutral dipole pairs, whose density and size, in turn, decrease further as T increases. We 
note that this picture is confirmed by our simulations — see the plot of the vortex density as a 
function of temperature on Fig. 5 (vortices in the spin model are identified with VF-bosons). 

Many — but not all — aspects of this picture can be studied analytically in the case of SU (2) 
QCD(adj) via the renormalization group equations for the corresponding electric- magnetic 
Coulomb gas. In particular, various critical exponents can be exactly calculated for QCD(adj) 
with SU (2) gauge group, due to the special properties of the appropriate electric-magnetic 
Coulomb gas [1]. The gas turns out to be dual to a XY-spin model with a ^-preserving 
perturbation and the phase transition in the spin model is in the same universality class as 
the deconfinement transition in 577(2) QCD(adj). 

In the case of SU (3) QCD(adj), a Coulomb gas picture of the theory and a spin- model dual 
to the Coulomb gas also exist. However, the nature of the small-L deconfinement transition 
for SU (3) gauge group is not well understood, as opposed to the SU (2) case, as we now 
describe. The results from [1] provide some guidance. An important property of the SU (3) 
QCD(adj) theory 5 on M 2 x x is that it exhibits electric-magnetic duality (valid in the 
temperature range (1.2)). The duality interchanges electric and magnetic fugacities as well as 
the electric and magnetic coupling (re or re -1 ) and takes re — > ^ [1]. By eqn. (1.5) of Section 
1.3.2, this is also a high-T/low-T duality, suggesting that if there is a single transition, the 
self-dual point is the critical point. At present, we do not know if at the self-dual point (1.4) is 
described by a conformal field theory (this would be the case if the transition was continuous 
and is exactly what happens in the calculable case of SU(2)). In [1], it was shown that the 
leading order (in fugacities) renormalization group equation for re do, indeed, exhibit a fixed 
point at the self-dual point of the electric-magnetic duality. However, both 14^-boson and 
bion fugacities are relevant at the self-dual point — indicating that a strong-coupling analysis 
is necessary, even for small fugacities at the UV-cutoff. Thus, determining the order of the 
transition calls for strong coupling methods — in particular, the present work. 

In Sections 1.3.2 and 1.3.3, we describe the duality map of SU(3)/Z^ QCD(adj) on 
t 2 x S[ x Sp to an "affine" XY-spin model. The discussion is meant to only introduce the 
reader to the duality and we only make what we consider the most relevant points. For more 
details, we refer to [1]. 

1.3.2 Spin model and symmetries in "natural" (root-lattice) coordinates 

The spin model dual to SU(3)/1*3 QCD(adj) is the theory of two coupled XY-spins. The 
most natural description of the model is in terms of two compact variables whose periodicity 
is determined by the SU (3) root lattice: 

e x = (el, e 2 x ) = e x + 2^1 = e x + 2™ 2 . (1.3) 

5 Electric-magnetic duality in the Coulomb gases relevant for QCD(adj) on I 2 x x Si does not hold for 
SU(N C > 4) due to the composite nature of the magnetic bions [1]. We will not study N c > 3 here. 
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Here are the simple roots 6 and 03, which appears in (1.4) below, is the affine root of 
577(3). In this Section, we use root-lattice coordinates, as they provide a natural description 
of the allowed electric and magnetic charges, the discrete symmetries, and the operators 
representing various nonperturbative objects in the gauge theory. The affine XY-model is 
defined by a lattice partition function with: 

N c =3 N c =3 

-pH = T~ cos2Vi • - 0*) +EE ycos2(ai - St-i) ■ 9 X . (1.4) 

x;p,=l,2 i=l x i=l 

where V{ (i = 1,2,3) are the weights of the defining representation, explicitly defined in 
Footnote 6, x denotes the sites on a 2d square lattice, and jl — the two unit lattice vectors. 

The first term in (1.4) is essentially the model used by D. Nelson [15] to describe the 
melting of a two-dimensional crystal with a triangular lattice. There, the two-vector 9 X 
parametrizes fluctuations of the positions of atoms in a 2d crystal around equilibrium (the 
"distortion field") and k is proportional to the Lame coefficients. A vortex of the compact 
vector field 9 X describes a dislocation in the 2d crystal. The Burger's vector of the dislocation 
is the winding number of the vortex, now also a two-component vector. For a description 
of the order parameters and phase diagram, see [15]. We only mention that the melting of 
the 2d crystal occurs due to the proliferation of dislocations (vortices) at high temperature, 
which destroys the algebraic long-range translational order. 

Now, we turn to the 5£7(3)/Z3 QCD(adj) meaning of the spin model. 

The "spin waves" in (1.4) — the fluctuations of 9 X — are the duals of the two massless 
photons. The dual photons are sourced by magnetically charged excitations — the magnetic 
bions, whose image in the spin model is the potential term in (1.4), discussed below. 

The vortices in the spin model (1.4) are dual to magnetic charges and thus describe the 
electric excitations in the gauge theory — the W-bosons of mass ~ j^, which can be excited 
at T > 0, as described in the previous Section. There are three kinds of vortices of equal 
lowest action that contribute to the dynamics. These are described in [1] — and, much earlier, 
in [15] — and correspond to the three lightest degenerate W-bosons in 577(3) QCD(adj) at 
small-L. The spin-spin (as well as vortex-vortex) coupling k, which determines the kinetic 
term of the dual photons, is expressed via the four-dimensional gauge coupling g±{L), the size 
of the spatial circle L, and the temperature T: 

In the melting applications of the "affine" XY-model, there is an exact £7(1) x £7(1) 
global symmetry 9 X -^9 X + c. The £7(1) 2 symmetry forbids terms in the action which are 
not periodic functions of the difference operator. In the gauge theory, this is the topological 
£7(1) 2 symmetry of the two free dual photons. This symmetry is explicitly broken by the 



e For future use, we give the weights of the defining representation: v\ = (|, j^g), vi = (— §, 5^75)1 ^3 = 
(0, — ^g). The roots are differences of weights and are given by c?i = V\ — V 2 = (1,0), S.2 = v 2 — v-.i = (— | , 2)1 
and affine root 0-3 = U3 — v\ = (— |, — ^). The weights obey Vi ■ Vj = — | for i,j = 1, 2, 3. 
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fact that the gauge theory possesses magnetic monopole-instanton solutions (essentially 't 
Hooft-Polyakov monopoles or "twists" [13] thereof). 

In the spin model dual to QCD(adj), the U(l) 2 symmetry is explicitly broken, by the 
potential "external field" term in (1.4), to a Z3 x Z3 symmetry, to be defined below. This 
breaking of U(l) 2 is crucial (and is the main difference between (1.4) and the models describing 
melting) as it captures the effect of the magnetically charged excitations in the thermal gauge 
theory, the magnetic bions responsible for confinement. The strength of the "external field" 
in (1.4), y, is proportional to the bion fugacity in the gauge theory. The bion fugacity, 
see [1, 12, 14], is small at small L, and is given by (up to a numerical constant): 

^^ZV^7 e " S °' (L6) 



where So ~ <f(L) * s ^ e mon ac ti° n - I n SU(S), there are three kinds of magnetic bions of 
different magnetic charges (~<3j — &j, i ^ j) but equal fugacities. The Coulomb gas of these 
three kinds of bions is represented in the spin model by the three equal-strength "external 
field" terms in (1.4). Thus, the model has an additional symmetry, permuting the arguments 
of the three cosines, simultaneously in the kinetic and potential terms. This symmetry is 
inherited from the unbroken Weyl group in the center-symmetric vacuum responsible for 
the SU(3) — > U(l) 2 breaking. This symmetry remains unbroken at zero temperature and 
we do not expect it to break at T > (indeed, the results of our simulations confirm this 
expectation, as the two magnetizations m 1 and m 2 (1.17) and their susceptibilities behave 
identically, but would not be expected to if the exchange symmetry between different bions 
and W bosons was broken). 

As described in Section 1.3.1, in the temperature range (1.2) the gauge theory has 
an electric-magnetic Coulomb gas description as a gas of magnetic bions and VF-bosons — 
both come in three varieties — interacting via electric and magnetic Coulomb potentials and 
Aharonov-Bohm interactions. A Coulomb gas description also holds for the spin model and 
is the reason behind the duality; see [8] for a lattice derivation and [1] for a continuum 
description. 

The two Z3 symmetries 7 of the spin model (1.4) are mapped to the topological Z 3 of 
SU(3)/li3 QCD(adj), associated with the nontrivial ~ki{SU{?>)/'L^), and the Z3* discrete sub- 
group of the anomaly-free global chiral symmetry. Note that while the nj-adjoint flavor 
theory has a larger chiral symmetry group, only Zg X is relevant for the deconfinement tran- 
sition on M 1 ' 2 x S^; the continuous SU{nj) chiral symmetry remains unbroken at small L. 



7 The reader not familiar with the Lie-algebraic constructions used in this Section is advised to proceed to 
Section 1.3.3, where a description in terms of coordinates rectifying the root lattice is given. While the physical 
interpretation of the two Z3 symmetries is not obvious in the "simulation" coordinates of Section 1.3.3, the 
presence of two Z3 symmetries is evident. 
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The topological 8 Z 3 acts on 9 by shifts on the weight lattice: 9 

Z 3 : 6-+ 9- 2-rrkui, k = 0,1,2; for any % = 1,2,3 . (1.7) 
The action of the other symmetry, the discrete chiral Z 3 X , is: 

z dx. e^d+^p, 3 = 0,1,2, (1.8) 

where p is the Weyl vector 10 of SU (3). It is useful to describe the action of the two symmetries 
on appropriate order parameters. For the Z 3 symmetry (1.7), the exponentials e 2%8 ' Uj can be 
taken as order parameters: 

Z| : e ^ ^ e 2i§^ _ (L9) 

To obtain the result on the r.h.s., one uses the definitions from Footnote 6. An insertion 
of e x ' v i in the path integral of the spin-model (1.4) corresponds to the creation, at x, 
of a pointlike monopole of minimal charge allowed by Dirac quantization in an SU (3) /Z3 
gauge theory. Recall that in SU(N c )/7 i n c theories, dynamical 't Hooft-Polyakov monopoles 
have magnetic charges taking values in the root lattice. Thus, their charges are larger than 
the minimal ones allowed by Dirac quantization (the minimal allowed magnetic charges take 
values in the finer weight lattice). 11 As explained in [1, 12], the operator (1.9), when lifted to 
the gauge theory, creates a minimal charge 't Hooft loop at x £ M 2 , wrapping around S^. 

In the S£/(3)/Z3 gauge theory, no electric probes in the fundamental representation are 
allowed. Thus, there is no center symmetry and the symmetry 12 associated with confinement 
is the "dual center" Z 3 symmetry associated with m(SU(3) Confinement is thus probed 
by the behavior of correlators of minimal-charge monopoles. At low temperatures, the Z 3 
symmetry is broken and the correlator of two separated minimal-charge monopoles (1.9) (or 
't Hooft loops around Si) should approach unity at large separation. Conversely, in the 
deconfined phase, the correlator has an exponential fall off and approaches zero at large 
distances, owing to the confinement of minimal charge monopoles at high-T. 

Under the discrete chiral Z 3 X symmetry (1.8), the operators creating dynamical mag- 
netic monopole-instantons — 't Hooft-Polyakov monopoles with charges in the root lattice — 
transform as: 

Z 3 X : e 2i °- d > 1 e i? ?e 2 *^ , (1.10) 

where (1.8) and Footnotes 6 and 10 were used to obtain the r.h.s. (it is easily seen that the 
operators (1.10) are invariant under Z 3 ). Monopole-instantons are charged under the discrete 

8 In [12] this symmetry is called "dual center symmetry", while earlier both "topological global Zn symme- 
try" [16] and "magnetic Zjv symmetry" [17] have been used. 

9 Despite appearances, there is only one Z3 symmetry in (1.7) — the difference between transformations with 
different i is a shift of 9 by 2irk times a root vector, which is an identification, as per (1.3). 

10 The Weyl vector is defined as p = /ti + (I2, where fli are the fundamental weights implicitly defined via 
o.i • flj = J £-. In the basis of Footnote 6, p = (|, ^). 

J1 See Sect. 2.2 in [12] for a rather formal (but valid for general simple Lie groups) discussion in the R 1 ' 3 x §^ 
context. Alternatively, Sect. 2.2.1 in [1] contains a less Lie-algebra loaded explanation for SO(3) and SU(2). 

12 See [18] for discussions of "center symmetry breaking at high-T", Zjv bubbles, and spatial 't Hooft loops. 
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chiral symmetry due to the intertwining of the topological shift symmetries and discrete 
chiral symmetries. This intertwining occurs due to the presence of fermion zero-modes in 
a monopole-instanton background, required by the Nye-Singer index theorem [19]. Note 
that magnetic bions have no fermion zero modes and are chiral-symmetry neutral, hence the 
external field term in the spin model (1.4) representing bions is neutral under (1.8). At zero 
temperature, the Zg X symmetry is broken, hence, the correlator of two monopole operators 
(1.10) approaches unity at large separation. We expect this symmetry to be restored at high 
T. In fact, we will see that — as was also found [1] in the SU (2) theory — the restoration of the 
topological and discrete chiral symmetries occurs at the same temperature, as we see only a 
single transition. 

Finally, we note that operators representing electrically charged objects — the VF-bosons — 
can also be defined in the spin model as "vortex creation operators" [1]. We have not simulated 
correlators of such operators — both because they are not local in the 6 X basis (and thus require 
significant more computer resources to study) and because they are not charged under any 
global symmetry of the SU(3)/Z^ theory — and will not discuss them here. 

1.3.3 Spin model and symmetries in "simulation" coordinates 

In the simulation, we find it more convenient to use coordinates that rectify the root lattice. 
An equivalent description of (1.4) is given in terms of two compact scalar fields on a 2d square 
lattice, <j) l x = (f>\, + 2tt (i = 1,2), which are both periodic by 2ir. These are related to 6 X as 
follows: 

ti = -0l + ^=%, & = -0l-±%. (1.11) 
The action of the lattice model (1.4) becomes, using V ^x = 4>x+H — 4>x'- 

- (37i = i^ E ( cos ( V /^) + cos (V^ x ) + cos(V^ - v^D) 

x;A=l,2 (L12) 

+ V^2 ( c °s(3^) + cos(3^) + cos(3(0* - 4> 2 X ))) . 

x 

We now define a "temperature" T", by ^ = proportional to the gauge theory temperature, 
see (1.5), as well an "external field strength" h = see (1.6). The lattice spacing of the 
spin model is taken to be equal to unity (and, physically, is of order L, the inverse UV cutoff 
scale of the 2d Coulomb gas). Thus, we arrive at the Boltzmann factor which is actually 
evaluated at every iteration of the Metropolis algorithm employed in our simulations: 

e -PH = e -# 

H ^~ \ E (cos(V^)+cos(V^)+cos(V^-V^)) (1.13) 

x;#=l,2 

-hJ2 ( c °s(3<^.) + cos(3^) + cos 3(0i - <t>D) • 

x 

The action (1.13) has a x Z| symmetry, Zg acting solely on cf) 1 by a shift of times 
integer, and Zg, acting only on (f> 2 by a ^ shift. Using the map (1.11), it is easy to see that 
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the topological Z 3 maps to a linear combination of Z3 and Z|: 
nrt aX . j 1 1 ^vr/c 2 , 9 27r/c 



J 3 



+ — , ^f-— , A; = 0,1, 2, (1.14) 



while the chiral Zg X is identical to Z 2 : 



Zj*:* 1 -^ 1 , </> 2 ^ 2 + ^, ,-0.1.2. (1.15) 



The main goal in this paper is to study the symmetry realization as the temperature T 1 
and coupling h in (1.13) vary. Below, we define the quantities measured in our simulations. 
The "magnetization" , M 3 , for each field component (j = 1 , 2) is the average j-th "spin" , a 
planar unit vector represented as e*^, over the entire lattice: 

Af' = 5>^ , i = 1,2. (1.16) 

X 

The order parameter for the Z3 symmetry is the mean magnetization m 3 per spin: 

Here and in Section 2 we use iV to denote the lattice width (the number of lattice points is 
iV 2 ). We also define the magnetic susceptibility x J of each component: 

X J = U N2 ^) U , i = l,2, (1.18) 

as well as the specific heat: 

We also measure the average energy per site (or average energy per spin, remembering that 
the spin is a two-component vector), defined as: 

E = w- (L20) 



2. Monte Carlo simulations 



We have simulated (1.13) on a square lattice of widths N = 8, 16, and 32, and with per- 
turbations of strength h = 0.1 and 1.0 using the Metropolis algorithm. These values for N 
and h were used for most of the direct measurements. In our study of the finite-size scal- 
ing of the susceptibility and specific heat we also used h = 0.5 and added N = 24 for all 
values of h. The data are obtained after 20000 Monte Carlo sweeps (a sweep consists of iV 2 
Metropolis iterations), taking data every 10 sweeps and disregarding the first 2000 sweeps for 
equilibration. Data acquisition is performed either by starting from a random configuration 



- 10 - 



at high temperature and then decreasing the temperature, or by starting from an ordered 
configuration at low temperature and then increasing the temperature. We take the temper- 
ature difference AT' to be as small as 0.01 in the critical region, and AT' = 0.001 in the 
case of finite-size scaling. Data points and error bars were generated from raw data using the 
bootstrap method 13 using 1000 resamples. 

2.1 Direct measurements: magnetization, energy, and vortex density 

We begin by a discussion of the measurements of magnetizations, energy, and vortex density. 
The order parameter of the system (magnetization per spin, eq. (1.17)) is shown on Figure 
3 as a function of the temperature T' for fugacities h = 0.1 & 1.0, for lattices of width 
N = 8, 16 & 32. We found that the magnetizations m 1 and m? behave identically and 
only show results for m . In all cases, the magnetization at high temperatures (beyond the 
observed transition) depends on the width of the lattice, with larger lattices more closely 
approaching the continuum limit of m l — > as T' — > oo. The transition also appears to 
occur more sharply as ./V increases (for a given h), and as the strength of the perturbation h 
increases. 

Upon running simulations with increasing and decreasing temperature, we have found no 
evidence of hysteresis. While this observation alone would favor a continuous transition, such 
a hypothesis is not consistent with our other findings discussed below (perhaps, the failure 
to detect hysteresis is due to a combination of: using rather small volumes, a large tunneling 
rate, and a sufficiently long simulation time). 

The energy per spin (1.20) vs. temperature is shown on Figure 4 for the various strengths 
of h and lattice widths. Similarly, the transition becomes sharper as h increases. The observed 
low-temperature limits of E = \jt = —3.3 & —6 for h = 0.1 & 1.0, respectively, correspond 
exactly to the state where all spins are aligned, as expected from eqn. (1.13). 

The vortex number was determined by calculating the change of the angle over each 
plaquette, from one site to the next, making sure that the angle difference between neighboring 
sites lies in the range from — ir to it. The vortex density of each field was recorded after every 
sweep. Vortices were detected with a sensitivity of 10~ 5 , meaning that at each plaquette, 
see [21], a vortex of type 14 (1, 0) — i.e., a vortex of (ft 1 but not </> 2 — was detected if the sum of 
the difference in the angles (p 1 of the spins around the plaquette's border was within 10~ 5 of 
2-7T (higher-order vortices are assumed, and are, in fact, suppressed). 

Figure 5 shows the evolution of the (1,0) vortex density as a function of temperature 
for the various strengths of perturbation and lattice widths. In all cases, the vortex density 
transitions from at low temperatures to about 0.11 at temperatures above the critical 
temperature. Again, this transition is sharpest for h = 1.0. The behavior of the vortex (i.e., 
VF^-boson) density is consistent with the interpretation of the deconfinement transition as 
occurring due to the "liberation" , upon increase of T, of the electric charges confined in small 
dipole pairs at low temperatures. 

13 For succinct explanation of the bootstrap method, see [20]. 

14 The density of (1,0) vortices is identical to that of (0, 1) vortices; see [1] for an review of the vortices in 
the multi-spin XY-model. The third vortex of the same fugacity, the (1, 1) vortex is difficult to distinguish 
from overlapping (1,0) and (0, 1) vortices in the simulation. 
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2.2 Magnetic susceptibility and specific heat 

The magnetic susceptibility (1.18) as a function of temperature, for the same systems as 
discussed in Section 2.1, is shown on Figure 6. In all cases, the susceptibility is roughly at 
both high and low temperatures, and attains a peak around the critical region of temperature. 
The height of the peak depends on the size of the lattice. In our simulations, the systems of 
lattice width N = 32 reached the highest peaks in susceptibility. The location of these peaks 
provide an estimate of the critical temperature. 

The specific heat (1.19) of the systems in question (h = 0.1, 1; N = 8, 16,32) is shown 
on Figure 7 as a function of temperature T' . The specific heat transitions from C = 1 at low 
temperatures to a sharp peak in the critical temperature region and subsequently decreases 
at high temperatures. The height of the peak varies with lattice size and increases with the 
strength of the perturbation h. We caution the reader against using the plots of this Section 
to infer details of the finite-size scaling — a more accurate study employing much finer AT' 
intervals in the critical region, as well as an additional lattice width and coupling, is presented 
in the following Section. 

2.3 Finite-size scaling and energy probability distribution 

To perform finite-size scaling, we studied the behavior of the magnetic susceptibility and the 
heat capacity in the critical region of each lattice width. We used additional simulations with 
AT' = 0.001 in the critical region for each lattice size. In addition to the two cases h = 0.1 
and h = 1.0, we also considered h = 0.5. In summary, in our finite-size scaling study, four 
different lattice widths N = 8, 16, 24, 32 were used for all three values of h = 0.1, 0.5, 1. 

The best finite-size scaling fit we found is consistent with the hypothesis of a first-order 
deconfining phase transition. As explained in [22] and references therein, in a finite volume 
the delta-function like behavior of the susceptibility x characteristic for a first-order transition 
is smeared out to a peak scaling as iV 2 (in d = 2) of width ~ N~ 2 . Similarly, the maximal 
value of the heat capacity C increases like N 2 . The maximal value of \ % as a function of T", 
Xmax(N) is expected to behave as: 

xla X (N) = c N 2 + . . . , (2.1) 

where dots denote terms subleading at large iV. 15 The maximum value of the specific heat 
as a function of T' at each N, C max (N), obeys: 

C_(A0 = ^^iV 2 + ... , (2.2) 

where T' c {oo) is the infinite volume critical temperature [22]. Here, E+ and E- are the average 
energies per site of the symmetric (E+) and broken (E-) phase at T' c . In particular, E + — EL 
is the latent heat (per site) of the transition. As the variation of T' C (N) with N we find is 
slight (thermodynamic fluctuation theory predicts ^^i^f^ ~ N~ 2 ), we can use the slope 

15 If one uses the usual definitions of critical exponents in a first-order transition to characterize the finite-size 
scaling, the scaling in (2.1) would imply "1 = d. 
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of the C m ax vs. A^ 2 line in (2.2) (found in our relatively small volumes) to estimate the latent 
heat per spin. 

On Figures 8, 9, and 10, we show, for h = 1, h = 0.5, and h = 0.1, respectively, 
the scaling of the maxima, with respect to the temperature, of x l iX 2 , an d C with N for 
N = 8, 16, 24, 32. The values of Cmax, Xmax an d Xmax are determined by fitting the raw data 
for each quantity with a degree 4 polynomial weighted with the error bars of the data (thus 
higher error bars have lower weights) and taking the maximum of that polynomial. The error 
bars in C max ,x ma x an< ^ Xmax are the two-sigma (95% confidence level) errors of the curve 
fitting. Figures 8, 9, and 10 indicate that the data for all iV and h are fitted well with the 
linear dependence on the volume expected of a first order transition. 

As a further check on the conjectured first order nature of the transition, we studied the 
probability distribution of the energy per spin. The histogram we plot was generated in a 
Monte Carlo run using 5700 Monte Carlo sweeps (after waiting 4000 sweeps for equilibration 
and then measuring the energy after every sweep) for A" = 32, h = 0.5 for several values of T", 
of which we show only three. Two distinct peaks are clearly visible in the energy probability 
distribution for T 1 = T' c on Figure lib. Conversely, only one peak is present at either below 
or above T[,, with its energy shifted in the right direction, see Figures 11a and 11c. The order 
of magnitude of the distance between the peaks is consistent with the estimate of the latent 
heat obtained using the slope in (2.2): taking T' c ~ 0.957 for h = 0.5, we find, using (2.2) and 
Figure 9(c), E+ — E- ~ .38. This value is consistent with the distance between the peaks, 
about 0.4, estimated from Figure 11, and provides a check on our data and the first-order 
transition hypothesis. 

For h = 0.1, on the other hand, we could not detect distinct peaks in the energy distri- 
bution and we do not show these results. This lack of distinct peaks could be because the 
h = 0.1 phase transition is a very weak first-order one, with small latent heat, 16 or because 
for small h the transition is continuous with 7/V smaller but close to 2. The precision of our 
data does not allow us to be definite on this issue, which requires further work to settle. 

Finally, we can test another prediction of the thermodynamic theory of finite-volume 
scaling in first order transitions. As Ref. [22] showed, for sufficiently large volumes and 
sufficiently close to T c , the heat capacity N z is expected to be a universal function of 
(T' — T{.(N))N 2 — in other words, when plotted, the curves for different A^ should collapse. 
To test this, we plot C(N, T')N~ 2 as a function of (T" - T{.(N))N 2 for N = 8, 16, 24, 32, for 
h = 0.5, on Figure 12. The collapse of the curves for different A" (apart from the smallest 
N = 8) onto each other is reasonably good for T' < T' c , but both the statistics and the 
collapse are poorer in the high-temperature phase. 

In summary, we have presented ample evidence that for moderately large values of the 
bion fugacity (for sure, h > 0.5, but we have not simulated any h between 0.1 and 0.5) 
the deconfinement transition is a first-order one. Larger values of the bion fugacities are 
expected when the decompactification limit L — > 00 at fixed strong coupling scale A is ap- 
proached, as the bion action becomes smaller upon decompactification. Precisely a first-order 

16 Indeed, the slope on Figure 10(c) is much smaller compared to the one for h — 1 or 0.5 with similar T c . 
Assuming a first order transition and thus eq. (2.2), we can estimate the latent heat to be ~ 0.14. 
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deconfinement transition was found in lattice simulations of the full QCD(adj) nj = 2 gauge 
theory [23]. In the next Section we will summarize our findings, outline avenues for future 
studies, and speculate on the phase diagram in the /3-L plane for various values of nt. 



3. Conclusions and outlook 

3.1 Summary of results 

This paper is devoted to a study of the phase transition in the "affine" XY-spin model dual 
to the 5C/(3)/Z3 QCD(adj) theory on I 2 x §J x S^, at small L, with a particular emphasis 
on the nature of the deconfining phase transition. In our Monte Carlo study, we found that 
as the temperature is increased, there is a single deconfining and discrete-chiral restoration 
phase transition, which is discontinuous at least for h > 0.5, as found for the deconfinement 
transition in SU(3) 4d QCD(adj) [23] (with nj = 4). Finite-size scaling of the susceptibility 
and heat capacity and a study of the probability distribution for the energy were used to 
corroborate our conclusion. 

We now continue with comments on possible future work and speculations. 

3.2 Further studies of (selfdual) electric-magnetic Coulomb gases 

The partition function of QCD(adj) on M 2 xSi x§^ in the regime (1.2) can be cast as the grand 
canonical partition function of a multi-component Coulomb gas of electrically and magneti- 
cally charged particles interacting through logarithmic (dual) Coulomb and Aharonov-Bohm 
phase interactions. Schematically, it can be written as: 



z = £ S / ILj dxidyje 



'^bion S\V ^W/bion 



(3.1) 



Here, y and y are the VF-boson and bion fugacities and the sums/integrals are over arbitrary 
numbers and spatial distribution of the various kinds of VF-bosons and bions. The actions 
Sbiom Sw, an d S\Y/bi on contain their electric and magnetic Coulomb and Aharonov-Bohm 
interactions (an explicit expression for these can be found in Eq. (4.9, 4.11, 4.13) in [1]). 

We show Eq. (3.1) here only to note that the electric-magnetic Coulomb gas parti- 
tion function (3.1) can be given a different description, which may be useful for analyti- 
cal studies. To write it down, we introduce a two-compoenent, two-dimensional scalar field 
$ = $2) and its dual field = 62)- These obey the equal time commutation rela- 
tion [0*(x), = —i5 l:1 9(x — y), where 6 is the Heaviside theta function, so that II^ = d x Q 
is the momentum conjugate to see, e.g., [24]. Then, the 2d Coulomb gas partition function 
(3.1) can be represented as the Euclidean vacuum functional of a one-dimensional quantum 
system with Hamiltonian density: 



H 



~M) a + ~M) a 



N c =3 

£ 

i=l 



y cos 





9 



cti - di-x)$ 



+ y cos 



9_ 



.(3.2) 



Here, y and y are the magnetic and electric fugacities also appearing in (3.1), g = gi(L), 
L, and T are the gauge theory quantities appearing in (1.5,1.6). The deconfinement phase 
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transition in the SU(3)/1a3 QCD(adj) theory is the quantum phase transition in the theory 
with Hamiltonian density (3.2). As already advertised, the theory exhibits electric-magnetic 
duality, which, in the variables of (3.2), takes the simple form: 17 

^ T ~ Q Att\/3LT 
$ 2 *+ei, $i^9 2 , y^y, -7=^— • (3.3) 



The last relation in (3.3) is the already mentioned — symmetry, with k given in (1.5), 
and shows that this is also a high-T/low-T duality. At the self-dual point, it is easy to see 
that both the electric and magnetic potential couplings (i.e., the two cosines) in (3.3) are 
relevant (in stark contrast to the SU{2) QCD(adj) case [1], where they are marginal). 

Now that we have presented the selfdual sin-Gordon formulation (3.2) of the model, it is 
of interest to contrast the Monte Carlo results on the phases of the "affine" XY-spin model 
with previous analytic studies of electric-magnetic Coulomb gases in 2d. The leading order, 
in bion and H^-boson fugacities, renormalization group equations (RGEs) of the Coulomb 
gas dual to SU(3) QCD(adj) [1] have a fixed point for the electric (and magnetic) coupling 
k. The fixed point occurs at the self dual point = \/l2 ~ 3.46. However, this analysis 
was (rightly) not considered conclusive — at «*, the bion and V7-bion fugacities are both 
relevant, indicating that a higher-order (indeed, strong coupling) analysis is necessary. Our 
present study indicates that no fixed point is present, at least for sufficiently high fugacities; 
instead, one finds a first-order phase transition. These findings would also be in contrast with 
conclusions drawn by a straightforward extension of an earlier analysis of RGEs to the next 
nontrivial order in fugacities. In Ref. [25], a similar self-dual Coulomb gas associated with the 
SU(3) algebra was studied, and the RG analysis was performed using the selfdual sin-Gordon 
representation similar to (3.2). The model studied there is not identical to our: as opposed 
to (3.2), in [25] only c?j's appear in the first cosine, rather than a, — accounting for 

this difference is trivial and allows us to recast the results for the RGEs of [25] to our model. 
Taking them at face value, one might be tempted to conclude that there is a fixed point; 
however, the approximations used in [25] to derive the RGEs also do not strictly hold for 
(3.2). The lack of fixed point that we find illustrates the danger of relying on an uncontrolled 
approximation to infer the existence of a fixed point. 18 



Invariance of (3.2) under (3.3) follows after using the expressions for the ST/ (3) affine roots <5i given in 
Footnote 6. Note that a Hamiltonian identical to (3.2) describes also the Coulomb gas for SU(N C ) QCD(adj), 
the only difference being that there are N c — 1 fields (<E>, Q) and that the sum over the affine roots now extends 
from 1 to N c . It is easily seen that for N c > 3 there is no analogue of the electric-magnetic duality (3.2); 
the physical reasons are explained in [1]. One consequence of this lack of self-duality is that for N c > 3 the 
couplings (~ g and ~ 1/<j) inside the two cosine-potentials in (3.2) renormalize differently. 

18 As already mentioned, the apparent lack of a fixed point in SU(3) QCD(adj) is in contrast with the 
SU (2) case, where a continuous transition was found. A line of self-dual fixed points extends to weak coupling 
(small fugacities), thus allowing a controlled calculation of various critical exponents [1] (in fact, using 2d 
CFT methods, one can show that the phase transition in the spin model dual to SU(2) (adj) is second order 
for all values of the fugacities, see the first reference in [26]). The spin model dual to SU(2) (adj) case is 
a single-component XY-model with a Z4-preserving perturbation. Continuously varying critical exponents, 
corresponding to the fixed-line behavior, have also been seen in lattice simulations of the Z^-model [27]. 
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In the future, it would be of interest to further study (3.2), both for N c = 3 and higher, 
with an aim towards understanding the small-/i behavior (e.g., continuous vs. discontinuous 
transition) that we were inconclusive about. This could possibly be done by applying 2d CFT 
methods, as in [26], to (3.3) and is an interesting task for future work. 

There are also some obvious generalizations of the present Monte Carlo study. Notably, 
the SU(N C > 3) QCD(adj) theories on R 2 x S^j x are also dual to electric-magnetic 
Coulomb gases, given by the higher- N c generalizations of (3.2). Due to the "composite" 
nature of magnetic bions, these gases do not exhibit electric-magnetic duality for N c > 3 and, 
at present, we have no clues as of the nature of the deconfinement transition in these theories 
(the leading-order RGEs do not give any useful information and there is no self-dual point). 
It should be relatively straightforward, if more resource consuming, to pursue this question 
numerically. 

3.3 The global j3-L phase diagram 

We end with a discussion of the hardest and most interesting question: the relation of the 
small-L deconfinement transition to that in the large-L theory. We first recall that in SU (3) 
QCD(adj), "large-L" lattice studies have found that the deconfinement and continuous non- 
abelian and discrete chiral symmetry transitions do not coincide and that deconfinement is 
a first order transition [23]. 19 Upon increasing the temperature, the deconfinement transi- 
tion happens first, while chiral symmetry is restored only at higher temperatures. Based on 
this knowledge as well as our small-L studies, we have drawn, on Figure 1(a) and (b), two 
possible phase diagrams of SU(3) QCD(adj) — for rif such that the infinite-L theory is con- 
fining. Zero-temperature confinement occurs for rif < n c p (here nj~' is the critical number of 
Weyl adjoints above which the theory is conformal, presumably ~ 4 or 5 for QCD(adj)). For 
rif = 5, most likely a conformal case, the simplest possible phase diagram is given in Figure 2 
and the arguments leading to it are discussed further below. 

The phase diagrams on Figure 1 were drawn using the information that we have on the 
phases from large-L lattice studies and from our small-L semiclassical and numerical studies. 
For large-L, our phase diagrams assume distinct deconfinement and chiral transitions as 
in [23] and in drawing the phase diagrams we have assumed that only three phases can 
coexist (barring accidents) , since for exactly massless fermions there are only two parameters 
to vary, f3 and L. As already mentioned, the small-L dynamics of QCD(adj) always preserves 
the SU(rif) chiral symmetry. In the confining theories, the continuous (as well as discrete) 
chiral symmetry is believed to be broken in the large-L limit. Hence, an extension of the 
present methods to incorporate it is required — if one is to continue the small-L studies to 
larger- L — for generic rif, see [6] for discussion. 

There is one exception worth mentioning, however. For SU(N C ) with nj = 5 adjoint 
Weyl fermions, the theory is very likely conformal in the infinite-L limit due to the existence 
of a Banks-Zaks infrared fixed point at (reasonably — i.e., numerically, but not parametrically) 

19 Thus, the order of the small-L deconfinement transition for QCD(adj) is the same as the 4d theory 
transition. This is in contrast to deformed pure- Yang-Mills theory, where, in SU(N C ) theories it was found 
that the small-L transition is second order, in the universality class of Zn c parafermion 2d CFT, see [26]. 
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L L so 

(a) One simplest possible /3-L phase diagram for (b) The other simplest possible /3-L phase diagram 
n f < nj' . for n/ < nj\ 



Figure 1: The two simplest possible phase diagrams connecting the small- L and large- L behavior 
for theories which are confining in the infinite-!/ zero-T (infinite-/? limit). The thick red lines in the 
lower left hand corner denote the location of the discontinuous phase transition studied here. In both 
diagrams, here and in Figure 2, "C" and "D" refer to confined and deconfined phases, while "x" and 
"dx" denote the realization of the SU(rif ) continuous and Z 2 N c n f discrete chiral symmetries in the 
various regions. /3d and f3 x are the deconfmement and 5?7(n/)-restoration temperatures, respectively, 
in the infinite-!/ theory. L x denotes a critical radius L beyond which SU{nf) is broken (at zero 
temperature). 
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Figure 2: The n/ = 5 theory is (most likely) conformal with a, presumably, sufficiently weak-coupling 
Banks-Zaks like fixed point a*. If so, upon increasing L, the semiclassical analysis is always applicable. 
In particular, it shows that the conformal zero-temperature theory confines at any finite L with an 
exponentially small mass gap [6]. There is a thermal dcconfinemcnt and dx-brcaking transition with 
critical temperature given, for SU(2), by T C (L) — [1]. The transition is second order for SU(2). 
The line drawn on the figure corresponds to T C (L) drawn with the two-loop SU(2) coupling. The 
transition for SU(3) is — probably, see Section 2.3 — discontinuous and occurs at temperatures T' ~ 1, 

2 I T \ 

parametrically giving T C (L) - similar to the SU(2) case. The semiclassical theory, now expected 

to be valid at all L, predicts that the discrete chiral symmetry, but not SU(rif), also breaks at the 
deconfmement transition. 



weak coupling. In this case, the phase diagram considerably simplifies, as shown in Figure 2. 
For SU(2), the confined (and discrete chiral-broken) 20 and deconfined (and discrete chiral- 

20 We note that in the weak-coupling semiclassical regime, the order parameter for the "dx" ^2N c nf symmetry 
are the exponentials of the dual photons [3], e' a ' a , while at strong coupling, local 5'?7(n/)-singlet but "dx" 
nonsinglet order parameters can be built from multifermion operators like det A" a A^-, where a = 1 . . . N£ — 1, 

a, P — 1,2, and i,j = 1, are adjoint, SL(2, C)-Lorentz, and flavor indices, respectively. 
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symmetric) phases are separated by the line T C (L) = ^jp, where a = |-, up to exponentially 
small corrections, as was shown in [1]. The coupling a(L) is, to two-loop order: 



47T 

a(L) = - 
Po 



1 ^loglog(l/A 2 L 2 ) 



log(l/A 2 L 2 ) /3 2 (log(l/A 2 L 2 )r 



(3.4) 



where A is the strong coupling scale (in an IR conformal theory, this is the scale where 
the behavior of the running coupling switches from asymptotically free to conformal), (3q = 
(22 — 4nj)/3 and j3\ = (136 — 64ny)/3. Eq. (3.4), with nj = 5, was used to draw the 
phase transition line (3 C (L) = on Figure 2. We emphasize that the usual strong coupling 
problem mentioned in Section 1.2 is avoided because the coupling (3.4) runs, at large-L, into 
the fixed point at a* = — = — 0.136 rather than into the infrared "Landau pole" (the 
bion action at the fixed point is — ~ 46, and the semiclassical approximation is still expected 
to hold [6]). 

The global — i.e., valid at any L — phase diagram for QCD(adj) for nj = 5 on Figure 2 
can be viewed as a prediction of our semiclassical studies and could, in principle, be verified 
by lattice simulations. We end with the remark that, one, not surprising, conclusion is that 
the phase diagram of QCD(adj) is rather complicated and that, in this paper, we have only 
explored a small corner tractable by a combination of semiclassical and numerical methods. 
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Figure 3: Order parameter vs. temperature T' of the Z3 x Z3 model with symmetry-breaking 
fields of strength (a): h = 0.1, and (b): h = 1.0 in lattices of width N = 8 (green crosses), 
N = 16 (blue stars), and N = 32 (red circles). 
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Figure 4: Energy per spin vs. temperature T' of the Z3 x Z3 model with symmetry-breaking 
fields of strength (a): h = 0.1, and (b): h = 1.0 in lattices of width N = 8 (green crosses), 
N = 16 (blue stars), and N = 32 (red circles). 
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Figure 5: Vortex density vs. temperature T' of the Z3 x Z3 model with symmetry-breaking 
fields of strength (a): h = 0.1, and (b): h = 1.0 in lattices of width N = 8 (green crosses), 
N = 16 (blue stars), and N = 32 (red circles). 
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Magnetic susceptibility vs. temperature T" of the Z3 x Z3 model with symmetry- 
fields of strength (a) h = 0.1 and (b) /i = 1.0 in lattices of width N = 8 (green 
N = 16 (blue stars), and iV = 32 (red circles). 
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Figure 7: Specific heat vs. temperature T" of the Z3 x Z3 model with symmetry-breaking 
fields of strength (a) h = 0.1 and (b) h = 1.0 in lattices of width N = 8 (green crosses), 
N = 16 (blue stars), and iV = 32 (red circles). 
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(a) Probability distribution of energy for T < T' c . (b) Probability distribution of energy for T' = T' c . 
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(c) Probability distribution of energy for T' > T' c . 

Figure 11: Probability distribution of the energy per site for h = 0.5, iV = 32 for subcritical, 
critical and supercritical temperatures. The double peak distribution at T' c indicates phase 
coexistence and a first-order phase transition. The distance between the peaks is consistent 
with the latent heat calculated from the slope of the C ma x{N) vs. N 2 line, see the main text. 
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Figure 12: A plot of C(N, T')N~ 2 vs. (T" - T' C {N))N 2 . The data for different N shows a 
reasonable collapse, except for the smallest N = 8 volume. 
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